Competing Supersolid and Haldane Insulator phases in the extended one-dimensional 

bosonic Hubbard model 
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The Haldane Insulator is a gapped phase characterized by an exotic non-local order parameter. 
The parameter regimes at which it might exist, and how it competes with alternate types of order, 
such as supersolid order, are still incompletely understood. Using the Stochastic Green Function 
(SGF) quantum Monte Carlo (QMC) and the Density Matrix Renormalization Group (DMRG), we 
study numerically the ground state phase diagram of the one-dimensional bosonic Hubbard model 
(BHM) with contact and near neighbor repulsive interactions. We show that, depending on the 
ratio of the near neighbor to contact interactions, this model exhibits charge density waves (CDW), 
superfluid (SF), supersolid (SS) and the recently identified Haldane insulating (HI) phases. We 
show that the HI exists only at the tip of the unit filling CDW lobe and that there is a stable SS 
phase over a very wide range of parameters. 
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Since its introduction in 1989 [Ij, the bosonic 
Hubbard model (BHM) has attracted continued interest 
due to its very rich ground state phase diagram 
especially in the presence of longer range interactions. 
Direct experimental relevance was established with the 
reahzation of this model Hamiltonian, with tunable 
parameters [^, in systems of ultra-cold bosonic atoms 
loaded in optical lattices [3|. In its simplest form, a 
single boson species with only contact repulsion, the 
system exhibits two phases in the ground state [l|, a 
superfluid (SF) and an incompressible Mott insulator 
(MI) depending on the particle filling and the interaction 
strength. Extensive quantum Monte Carlo (QMC) 
simulations have established that, when longer range 
interactions are included, the supersolid (SS) phase 
can form for a wide range of parameters and lattice 
geometries in one, two and three dimensions |4|-|l6|. 
Typically, the SS phase is produced by doping a phase 
exhibiting long range charge density order (CDW). 
However, further QMC work has shown that, remarkably, 
for a wide range of parameters, the extended BHM 
with near neighbor interactions exhibits a SS phase even 
at commensurate fillings in two and three dimensions 



fillings of the system or is it a special property of the unit 



In addition, it was shown that the one-dimensional 
extended BHM with next near and/or near neighbor 
interactions admits another exotic phase at a filling of one 
particle per site; the Haldane insulator (HI) [19|, |20| . The 
HI is a gapped insulating phase characterized by a highly 
non-local order parameter like the Haldane phase |2ll.l22| 
in integer spin chain systems (see below) . This gives rise 
to several questions. Does the HI exist for other integer 



filling case? The SS phase found in one dimension [11 1 
was obtained by doping a CDW phase: Does this phase 
also exist for commensurate fillings in one dimension for 
parameter choices similar to those in two |17| and three 
dimensions [18[? If the SS phase exists for commensurate 
fillings, where is it situated in the phase diagram relative 
to the CDW, MI and HI phases? The phase diagram 
at unit filling for the BHM with contact (U) and near 
neig hbor (V) interactions was determined via QMC 
[la m and found to have SF, MI and CDW phases 
but no SS. Subsequently, the (/i,i) phase diagram of the 
extended BHM, for a fixed V/U ratio, was obtained using 
Density Matrix Renormalization Group (DMRG) j25| . 
but showed only evidence for MI, SF and CDW. More 
recent work, also based on the DMRG, has found [26[ 
no SS phase in the ([/, V) plane at unit filling but the 
question of other fillings was not addressed. Reference 
[2q also found a HI phase sandwiched between the MI 
and CDW phases which was not present in [23|, |24|. As 
wc shall see below, the HI phase was not found in the 
earlier work because the superfluid density in this phase 
vanishes very slowly with the system size and the largest 
sizes that could be accessed at the time were 64 sites. 

Theoretical studies of this system using bosonization 
have led to mixed results. The HI was obtained 
and characterized with bosonization [20[ but consensus 
is absent on whether the SS phase exists in this 
model. Even though older studies did not specifically 
mention it [27| or even argued that it did not 
exist [25[, more recent studies seem to demonstrate the 
presence of the SS phase [28(, even without nearest 



neighbor interaction [29|, for both commensurate and 
incommensurate filhngs. However, the precise nature 
of order and the decays of the relevant correlation 
functions are still far from settled. For instance, some 
studies predict that the single particle Green function 
decays exponentially in the SS phase while the density- 
density correlation function decays as a power [SOj ; others 
predict that both of these correlation functions decay 
as powers [28[. Finally, the universality class of the 
transition to the SS phase remains largely unexplored. 

In this Letter we answer some of these questions using 
the Stochastic Green Function (SGF) QMC algorithm 
[31! and the ALPS ^ DMRG code to obtain the phase 
diagram of the extended BHM in one dimension, 

i i 

+V^n,n,+i. (1) 

i 

The sum over i extends over the L sites of the lattice, 
periodic boundary conditions were used in the QMC 
and open conditions with the DMRG. The hopping 
parameter, t, is put equal to unity and sets the energy 
scale, a,; (aj) destroys (creates) a boson on site i, n.i = 
a\ai is the number operator on site i, U and V are 
the onsite and near neighbor interaction parameters. 
All results presented here were obtained at the fixed 
ratio V/U = 3/4 which favors CDW phases over MI at 
commensurate fillings when U is large. 

Several quantities are needed to characterize the phase 
diagram. The superfluid density is obtained with QMC 
using [33| 
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where W is the winding number of the boson world lines, 
d is the dimensionality and (3 the inverse temperature. 
The CDW order parameter is the structure factor, S{k), 
at fc = TT where 
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and the momentum distribution, n^, is given by 
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The charge gap is given by Ac(n) = /i(n) — ^(n — 1); 
the chemical potential is fJ-{n) = EqIu + 1) ~ EQ{n) 
where Eo(n) is the ground state energy of the system 
with n particles and is obtained both with QMC and 
DMRG. The neutral gap, A„, is obtained using DMRG 
by targeting the lowest excitation with the same number 
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FIG. 1: (color online) Shows several quantities for p = 1 
as functions of t/U with the fixed ratio V/U = 0.75. (a) The 
CDW order parameter for L — 6A and ps for L = 64, 100, 150; 
(b) the parity and string order parameters; (c) the neutral and 
charge gaps, (a) and (b) were obtained with QMC and (c) 
with DMRG. The region between the vertical dashed lines is 
the HI. In the HI, ps — ^ very slowly as L increases while the 
string parameter is essentially constant. Note the difference 
between the neutral and charge gaps. The gaps are given in 
units of the hopping i. 



of bosons. In both CDW and HI phases, the chemical 
potentials at both ends are set to (opposite) large enough 
values, in DMRG, such that the ground state de generacy 
and the low energy edge excitations are lifted [l9|, l25|. 
With the SGF we did simulations in both the canonical 
and grand canonical ensembles. 

For large values of U and V at p = I, a site typically 
has rir = 0, 1,2 particles with higher occupations being 
very rare. The s yste m then becomes analogous to an 
5 = 1 spin chain [l9[ with Sz{i) ^ Srii = ni — p [p = 1 
here) taking values of 0, ±1. Consequently, string and 
parity operators can be defined [191, |20J to characterize 
the Haldanc Insuating phase, 

a(N-.?|) = (<5n,e'^'S'--'"^fcj), (5) 

Op{\i-j\) = {e-'^L^'"'}, (6) 

where = tt for S = 1. The corresponding value of the 
order parameter is obtained in the limit |i — j| ^ 00; in 
practice we take the order parameters to be Os/p{Lmax) 
where, in QMC with PBC, L^^x = L/2 and in DMRG, 
with OBC, Lmax is the longest distance possible before 
edge effects start being felt. For higher integer filling, 
p = 2,3 . . ., Q ^ -K and has to be determined as discussed 
in [34|. It was shown |19l . l20| that the phases are 
characterized as follows [37[: In MI, ps = 0, ^(Tr) = 0, 

A, = A„ ^ 0, Op(i™ax) ^ 0, Os{Lraax) = 0; iu CDW, 
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FIG. 2: (color online) Same as Fig.[T]but at p = 2. (a) QMC 
simulations show that in the interval between the two vertical 
(black) dashed lines there is similtaneous SF and CDW order 
and, therefore, a supersolid phase (SS). The vertical (red) 
dot-dash line is where the L — >■ oo extrapolated neutral (A„) 
and charge (Ac) gaps vanish in DMRG (c). The CDW-SS 
transition is between t/U = 0.33 (QMC) and t/U = 0.355 
(DMRG). The difference between the two values could be due 
to the difference in the boundary conditions, open for DMRG 
and periodic for QMC. (c) also shows the L ^ oo extrapolated 
CDW order parameter, right (red) triangles, and the Fourier 
transform of {ni){nj), left (black) triangles, obtained with 
DMRG to probe the disappearance of CDW order. Both 
DMRG and QMC give the SS-SF transition at t/U « 0.425. 
Note that, unlike Fig. [T] the charge and neutral gaps (c) are 
essentially always the same. 



Ps = 0, Sin) ^ 0, Ac and A„ ^ 0, 0.p{Lmax) and 
a(ima.) ^ 0; in SF p, ^ 0, S{ix) = 0, A,. = A„ = 0, 
O^iU^a:,) -= Os{L„,ax) = 0; in HI, p, = 0, S{tt) = 0, 
Ac and A„ ^ 0, Op{L„^ax) = 0, 0,{Lmax) ^ 0. In the 
SS phase, CDW and SF order coexist and one expects 
Ps ^ 0, S{7t) ^ 0, Ac = A„ = and, because ^(Tr) =^ 0, 
Op{Lmax) ^ 0, Os{L,nax) ^ 0. The gaps, Ac and A„ 
behave in a subtle way in the CDW and HI phases (see 
below). 

We start with the system at p = 1 and study the phases 
as functions of t/U. Figure [1] shows the dependence 
on t/U of Ps, S{k — 7r), Os, Op, Ac and A„. Figure 
[TJc) shows that in the CDW phase Ac > A„ and that 
A„ = 0, Ac 7^ at the CDW-HI transition. The HTSF 
transition is signaled by Ac = A„ — >• [19|, l20| . Finite 
size scaling of the DMRG results show that A„ -> at 
t/U « 0.32 ± 0.01. Therefore, according to the criteria 
discussed above, the system is in the CDW phase for 
t/U < 0.22 and in the SF phase for t/U > 0.32. For 
0.22 < t/U < 0.32 (between the two vertical dashed 
lines), the system is in the HI phase with ps —>■ as the 



size increases. Note how slowly ps ^ with increasing L 
and how insensitive Os is to the finite size. This makes 
Os a more reliable indicator at moderate system sizes. 
Our CDW-HI transition at t/U = 0.22 agrees very well 
with Fig. 1 in [26[. However, the value we obtain for the 
HI-SF transition, t/U « 0.32 does not agree with the 
schematic dashed line in that figure. 

As mentioned above, the behavior at p = 1 may be 
understood by making the analogy with S = I spin 
chains. The question arises then as to whether such 
an analogy between this extended BHM at p = 2,3... 
and S" = 2, 3 . . . spin chains is valid and also leads to HI 
phases. This question is addressed in Fig. [2] for p = 2 
which shows qualitatively different behavior compared 
to Fig. [Ij While for low t/U both cases exhibit CDW 
phases, the behavior of Ac and A„, calculated with 
DMRG, is strikingly different as seen in Fig. [2][c): For 
p = 2, Ac = A„ and finite size scaling shows that 
they vanish together at t/U « 0.36, which indicates 
that there is no HI in this case. This is consistent 
with the absence of the Haldane phase in S" = 2 spin 
chains l35|- Nonetheless, for this fUhng, the system docs 



exhibit another salient feature: Indeed, figure ^a.) and 
(c) show from both QMC and DMRG that when the 
gaps vanish, ^(Tr) remains non-zero while ps also takes 
a non-zero value. <S'(7r) and ps both remain non-zero 
for 0.33 < t/U < 0.425 indicating the presence of a 
supersolid phase. The CDW-SS transition is estimated 
to be at t/U « 0.33 from QMC and t/U « 0.36 from 
DMRG while both DMRG and QMC give t/U « 0.425 
for the SS-SF transition. In Fig. U we show n^/L and 
S{k) in the SS phase at t/U = 0.35 and L = 64, 100, 128. 
We see that while rik/L -^ (sec center peak. A: = 0) as 
expected (since there is no condensate in one dimension) , 
the peaks in S{k) do not depend on L, indicating long 
range CDW order. This behavior is also confirmed by a 
finite-size scaling analysis of the DMRG results for sizes 
L = 64, 96, 128, 160. For the three phases (CDW, SS and 
SF), the CDW order parameter ^(Tr) is found to scale as 
S'o+S'i/i+S'2/i^, whereas uq/L is always found to decay 
as a power law ni/L". More precisely, one obtains the 
following results: for t/U ^ 0.25 (CDW): So = 3.4 and 
a « 0.95; for t/U = 0.39 (SS): 5*0 = 1.6 and a w 0.22; 
for t/U = 0.49 (SF): Sq « and a « 0.14. Note that, 
in both SS and SF phases, the parameter a is less than 
0.25, in agreement with a Luttingcr liquid description of 



the system [19|, |20|, |27|, l2g|. This scaling law and the 



insensitivity of ps to L, Fig. [UJa), confirm that this is 
indeed the SS phase. This surprising appearance of the 
SS phase at commensurate filling has also been observed 
in two and three dimensions [171 Il8l | . Furthermore, we 
find that the Green function, G{r) = (ajao), decays as a 
power in the SS phase with exponent « 0.5 at t/U = 0.34. 
For the present value of V/U, this behavior at p = 2 
is repeated at p = 3 (and presumably at higher integer 
fillings): As t/U is increased, the system goes from CDW 



1 

0.8 
^0.6 
~c 0.4 



m 2 



1- 





p=2! U=2.'857t, V/U=6.75 
Supersolid 



o-oL=64 
B-oL=100 




(b) 



-1 u "^ 

k 



FIG. 3: (color online) The dependence of the momentum 
distribution, rik/L and the structure factor, 5(fc) on the 
system size. rik/L ^ Q with increasing L while S{k) remains 
constant indiating long range CDW order. 
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FIG. 4: (color online) Main panel: The phase diagram of 
the BHM obtained with QMC and DMRG simulations. All 
results are for a system with V/U = 0.75. All symbols are 
from QMC and are for L = 64 sites, /3 = 128 (stars are 
for L = 128). Black lines near the tips of the CDW lobes 
are DMRC results for L — 192. Inset: A zoom of the tip 
of the p = 1 lobe. To the right of the vertical dashed line 
{t/U = 0.22) is the HI phase, to the left is the CDW. 



to SS to SF without exhibiting any HI phases. It appears, 
therefore, that, at least for V/U = 3/4, the analogy 
between integer spin chains and this extended BHM at 
integer fillings applies only at p = 1. 

The phase diagram in the {t/U, ji/U) plane for V/U = 



0.75 is mapped by calculating the charge gaps at 
commensurate fillings, multiples of L/2, and by making 
plots like Figs, m [5] which allow us to identify the various 
phases. We have also verified that, for incommensurate 
fillings, DMRG agrees with QMC on the boundary 
between the SS-SF. The result is shown in Fig.^l Figure 
|4] was obtained using QMC (all symbols) and DMRG 
(black lines near lobe tips). The QMC results are for 
L = 128 (stars) and L = 64 (aU other symbols), /3 = 128. 
The solid black lines near the lobe tips are obtained from 
DMRG with L = 192. The end points of the lobes are 
obtained by studying the finite size dependence of A„ 
using DMRG. The inset is a zoom of the tip of the p = 1 
lobe. 

Several comments arc in order. The p ~ 1/2 lobe 
is surrounded almost entirely by SF except for a small 
region of SS squeezed between it and the p = 1 lobe. The 
fact that in the extended BHM a SS does not exist when 
the p — 1/2 CDW phase is doped with holes, but does 
when it is doped with particles, was already addressed in 
111 ] . The p = 1 lobe sticks out of the SS phase and the 



part sticking out is, in fact, the HI phase. No other CDW 
lobe behaves this way. The p = 3/2 lobe terminates 
right at the boundary with the SF phase: To within 
the resolution of our simulations, the transition from 
the p = 3/2 CDW lobe goes directly into the SF phase 
without passing through the SS phase. This peculiar 
behavior for p = 3/2 was also observed with additional 
DMRG results for different values of V/U ranging from 
0.65 to 1: The SS layer between the CDW and SF phases, 
if present, is too thin to observe for the considered system 
sizes. An accurate determination of the (t/, V) phase 
diagram for this filling will require a more thorough finite 
size scaling analysis. All other CDW lobes, p > 2, are 
surrounded entirely by the SS phase. It is interesting 
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to compare this figure with Fig. 3 of 
mean- field predictions [36| . 

In this letter we examined the phase diagram of the 
extended BHM at a fixed ratio of the interaction terms, 
V/U = 3/4. Contrary to expectation, we found that 
this model at integer fillings does not always behave 
analogously to integer spin chains. In particular, only for 
p = 1 and at small t/U docs this happen and the system 
exhibits CDW, HI and SF phases. In the CDW phase 
at this filling, Ac > A„. At all other integer fillings, 
we found the HI phase to be absent and in its place 
a supersolid phase which indicates that the system at 
these fillings may not behave like an integer spin chain. 
Furthermore, for all CDW phases, except the one at 
p = 1, we found that A„ = Ac and that, unlike the p = 1 
case, both gaps vanish together as the CDW phase gives 
way to SS or SF. It is possible that, for a different V/U 
ratio, the SS-SF boundary will shift and cut the p = 3 
lobe (as it does in Fig. |4]with the p = 1 lobe) resulting in 
a HI phase. If this happens, it could mean there are two 
types of p = 3 CDW phases, one in which the neutral 



and charge gaps arc always the same (what we find here) 
and another CDW phase in which Ac > A„ as is the case 
for the p = 1 CDW. We have also shown that the single 
particle Green function decays as a power law in the SS 
phase. Finally, from a theoretical point of view, it would 
be interesting to characterize the universality classes of 
the different quantum phase transitions (CDW-SS-SF) 
occuring in the system. In addition, in order to have clear 
experimental signatures of the phases, in particular with 
cold atom gases which allow time-resolved measurements 
one would need to study the excitations of the system. 
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